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The process J/ip — * PP is studied using 57.7xl0 6 J/t/j events collected with the BESII detector at 
the Beijing Electron Positron Collider(BEPC). The branching ratio is determined to be Br(J/ij) — > 
pp) = (2.26 ± 0.01 ± 0.14) x 10" 3 , and the angular distribution is well described by d ™ g = 

1 + a cos 2 dp with a = 0.676 ± 0.036 ± 0.042, where 8 P is the angle between proton and beam 
direction. The obtained a value is in good agreement with predictions of first-order QCD. 
PACS: 13.25.Gv, 14.40.Gx, 13.40.Hq 
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I. INTRODUCTION 

The J/ijj meson is interpreted as a bound state of a 
charmed quark and charmed antiquark (cc). The decay 
process J/4 1 — * PP is an octet-baryon-pair decay mode 
that has been measured by DM2 DASP jj, Mark 
I 0, Mark II and Mark III @. 

In general, the angular distribution of J /ip — > BB can 
be written as 



dN 



dcos 9b 



oc 1 + a cos 6b, 



where 9b is the angle between the direction of the omit- 
ting baryon and the e~e + beam direction. Different the- 
oretical models based on first-order QCD give different 
predictions for the value of a [ESQ- Brodsky and 
Lepage assumed that the reaction is a process in the 
asymptotic region where all quark and baryon masses can 
be neglected. Claudson, Glashow, and Wise |7J included 
the baryon masses in their calculation, which significantly 
changed the previous prediction. Carimalo || consid- 
ered the effect of non-zero quark masses. The branching 
ratio and corresponding angular distribution will pro- 
vide some insight into the details of the baryon struc- 
ture H H, U 113 Chernyak and Zhitnitsky @ have 
investigated QCD sum rules for nucleon wave functions 
and calculated the branching ratio of the discussed de cay 
in terms of the wave function. Ping, Chiang and Zou [l(j 
studied the branching ratio as well as the angular distri- 
bution for J/ij) — * BB decay processes in naive quark 
model. Table [I] shows the measured branching ratio and 
a value of J/ip — > pp by previous experiments. Table ITTI 
shows the different theoretical predictions on a value. 



a value. The measurements are based on an analysis of 
57.7 xIO 6 ,//^ [prevents registered in the BESII detector 
at BEPC. 

BESII is a large solid-angle magnetic spectrometer 
that is described in Ref. [Tj. Charged particle mo- 
menta are determined with a resolution of a p /p=1.78% 



in a 40-layer cylindrical main drift cham- 

ber(MDC). Particle identification is accomplished by spe- 
cific ionization (dE/dX) measurements in the drift cham- 
ber and time-of-flight (TOF) measurements in a barrel- 
like array of 48 scintillation counters. The dE/dX reso- 
lution is o^e /,ix — 8.0%; the TOF resolution for Bhabha 
events is cttof — 180 ps. These two systems inde- 
pendently provide more than 3tr separation of protons 
from any other charged particle species for the momen- 
tum range relevant to J/ip — > pp decays. Radially out- 
side of the TOF counters is a 12-radiation- length bar- 
rel shower counter (BSC) comprised of gas proportional 
tubes interleaved with lead sheets. The BSC measures 
the energies and direction s of photons with resolution of 
<j e I E ~ 21%/v A B{GeV), ct = 7.9 mrad, and a z = 2.3 
cm. The iron flux return of the magnet is instrumented 
with three double layers of counters that are used to iden- 
tify muon. 

In this analysis, a GEANT3 based Monte Carlo pack- 
age(SIMBES) with detailed consideration of the detector 
performance (such as dead electronics channels) is used. 
The consistency between data and Monte Carlo has been 
carefully checked in many high purity physics channels, 
and the agreement is reasonable. 



II. EVENT SELECTION 



TABLE I: Previous experimental results for J/ip — > pp 



measured value about Br. and a 



Coll. 


jyobs 


Br(xltT 3 ) 


a 


Markl 


331 


2.2L0.2 


1.45±0.56 


Mark2 


1420 


2.16±0.07±0.15 


0.61±0.23 


Mark3 




1.91±0.03±0.16 


0.58±0.14 


DASP 


133 


2.5L0.4 


1.70±1.70 


DM2 




1.91±0.04±0.30 


0.62±0.11 


PDGfl4] 




2.12±0.10 


0.63±0.08(avg.) 



TABLE II: Different theoretical predictions on a value 



Theoretical values of a 


a 


authors 


1.0[6] 


Brodsky and Lepage 


0.4617] 


Claudson, Glashow and Wise 


0.69,0.70f8] 


Carimalo 



In this paper, we report high precision measurements 
of the branching ratio of J /ip — > pp and corresponding 



Events are selected with two and only two well re- 
constructed and oppositely charged tracks. The polar 
angle of the track, 9, is required to be in the range 
cos 6*| < 0.8. In order to remove the cosmic rays, the 
difference of the flying time between the positive and 
negative track, \t+ — t-\, is required to be less than 4.0 
ns. Information from the TOF system is used to identify 
protons. Both tracks are required to be not identified 
as a 7T or a K by the TOF, i.e. the measured flying 
time and the expected flying time assuming the track 
to be proton, pion, or kaon should satisfy the following 
requirement, \t meas - t exp (p)\ < \t meas - t exp (ir)\ and 
\tmeas - t exp {p)\ < \t me as - t exp (K)\. Meanwhile, since 
J r /tp — > pp is a two-body decay, the pp pairs are back- 
to-back; we require the two-track opening angle to be 
larger than 175°. Finally, for both tracks, the absolute 
difference between the measured momentum and its ex- 
pected value, 1.232 GeV/c, is required to be less than 
110 MeV/c. 

Figures fTlaVfd) show the proton momentum distri- 
bution for events surviving different selection criteria: 
Fig. nja) is after the particle ID and cosmic rays veto; 
Fig.^b) is after the back-to-back requirement; Fig.^c) 
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is the antiproton momentum after the same cuts; and 
Fig. ^d) is the proton momentum distribution after se- 
lecting the momentum of antiproton with ±110 MeV/c 
cut. From this last figure it is apparent that most of the 
background has been rejected. 



III. BACKGROUND ESTIMATION 

Two methods are used to estimate the background 
level. Method one uses a J/^ ^ anything MC sam- 
ple that is produced using the Lund-charm |l2j genera- 
tor. From this sample the level of background events in 
the data is determined to be about 1.5% of all selected 
events. This background is mainly due to J/ip — > jpp, 
J/ip — > ir°pp and J/tp — > jn c , n c — > pp. Before appli- 
cation of the proton momentum requirement, the back- 
ground level is about 1.7%. This study indicates that 
there is little background for high momentum values, but 
some at lower momenta, as can be seen in Fig.E](d). 

In method two, possible two-body background chan- 
nels, such as J/tp — » e + e~, J/ip — > J/ip — > 
K + K~, J/ip —* 7r + 7r~, J/ip — > ^pp, J/ip —> tt°pp, and 
J/ip — > jn c , n c —>■ pp are considered and generated ac- 
cording to the branching ratios listed in PDG(2002) [l4l |. 
The total background level for these events surviving the 
J/ip — > pp event selection criteria is about 0.7%. 

Since the background level is very small, we do not 
subtract the background but use 1.5% as the uncertainty 
of the background in the branching ratio measurement. 



IV. EFFICIENCY CORRECTION 

Since there are imperfections in Monte Carlo simula- 
tion, particularly when it simulates the tail of momentum 
resolution and TOF resolution, which will affect the effi- 
ciencies of momentum cuts and particle ID. In order to 
reduce systematic error as much as possible, the correc- 
tion on MC efficiency is needed. Since both momentum 
resolution and TOF resolution are cos 9 dependent, when 
we put cuts on momentum and TOF, the efficiency will 
depend on the angle 9. A reweighting method is used to 
do the efficiency correction. We define a cos 9 dependent 
weight factor (or correction factor) wtj (cos 9) as: 

wt J {cos9) = ef ta (cos9)/eJ lc (cos9), 

where j represents the jth selection requirement, 9 is the 
polar angle of p or p, and e stands for efficiency. When 
the weight functions are obtained, phase space MC gen- 
erator is used and 500,000 events are generated. Since 
the weight is obtained in each small bin of cos 9, the 
weight function will not depend very much on the actual 
angular distribution of generated events. To determine 
the weight for a given variable, the other, uncorrelated 
requirements are made more stringent to provide a pure 
sample. When the correlation among different selection 



criteria is so small that it can be ignored, the efficiency 
of the data can be expressed by the reweighted MC effi- 
ciency as: 

n 

e da ta(cos9) = e x H(e^ ta (cos9)) 

3=1 

n 

= e x JJ(ey* c (cos0) x wt j (cos9)) 

3=1 

n n 

= ex n s mc ( c ° s °) n wt 3 ( c ° s ^ 

3=1 3=1 

= e mc (cos0) x wt to t(cos9), 

where ex is the part of the efficiency which is not cor- 
rected; it is due to track reconstruction, geometry accep- 
tance, etc, and its uncertainty can be determined using 
another method(see Sec. IVllI below). The total correc- 
tion factor wttot is: 

wt to t(cos9) = 

wt p (cos9) x wtp(cos9) x wt P id(cos9) x wt pil i(cos 9), 

where wt p and wtp are the weights for the proton and an- 
tiproton momentum requirements, while wt P id and vutpid 
are the weights for the particle identification require- 
ments on proton and antiproton. 

V. ANGULAR DISTRIBUTION 

With e mc (cos6>) denoting to the efficiency obtained 
from Monte Carlo and wtt t(cos 9) for the total correction 
function of the efficiency, we fit the angular distribution 
of the data with the function f(cos9), 

/(cos 9) ~ A x (1 + acos 2 9) x e mc (cos6>) x wi tot (cos 9) 

where A is a constant. 

Fig. Ola) shows the angular distributions for 500,000 
phase space MC events and Fig. |2Jc) for the sideband 
background which is from p,p momentum distribution 
from 1.00 GeV/c to 1.10 GeV/c. Fig. E^b) shows the 
total weight curve and Fig.|2Jd) shows the result of fitting 
to the angular distribution of the data. The fit to the 
angular distribution of data gives 

a = 0.676 ± 0.036, 

where the error is statistical only. 

VI. BRANCHING RATIO OF J/ip -> pp 

After the final selection, the number of signal events is 
n°l s ta = 63316 

The Monte Carlo generator, which includes the angu- 
lar distribution of , dN a = 1 + a cos 2 9 V with a = 0.68 

a cos 9 P P 
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FIG. 1: The proton momentum distribution for events surviving different selection criteria: (a) after particle ID and cosmic 
ray veto; (b) after the back-to-back requirement; (c) the corresponding antiproton momentum distribution; (d) the proton 
momentum distribution after selecting the antiproton momentum as shown in(c). 



obtained from above fit, is used and 50,000 events are 
generated to estimate the selection efficiency for branch- 
ing ratio calculation. The MC-dctcrmincd selection effi- 
ciency e mc , including the effect of reweighting, is e mc — 
(48.53 ± 0.31)% and the branching ratio is determined to 
be: 

n obs 

Br(J/ip -> pp) = tot data = (2.26 ± 0.01) x 10~ 3 , 

N f/ip x e ™ c 

where Nj"^ is the total number of J/ ip in the data sample 
(57.7 x 10 6 ) [3. The error is statistical only. 

VII. SYSTEMATIC ERROR 

A. Systematic error of angular distribution 

When the fit parameters of the weight curve are 
changed by ler, the fitted value for a changes by 2.6%. 
This is taken as a systematic error. Another system- 
atic error from tracking reconstrucion is determined by 
using different MDC wire resolution models in the MC 
simulations. This affects the a value by 5.2%. From 
a p,p momentum sideband study, the uncertainty from 
background, including background level and background 
shape, is estimated to be 2.2%. The uncertainties 



from other sources such as from the fitting function of 
emc(cos#) are negligible. Adding these contributions in 
quadrature gives a total effect about 6.2%, i.e. 

a = 0.676 ±0.036 ±0.042, 

where the first error is statistical and the second system- 
atic. 



B. Systematic error of the branching ratio 

Systematic errors on the branching ratio measurement 
mainly come from MC statistics, the efficiency uncertain- 
ties of the particle ID, the momentum selection, the MDC 
wire resolution, the uncertainty in a, uncertainties in the 
background level, and the total number of J/tp events. 

The MC statistics gives a systematic error about 0.6%. 
As discussed in Sec. IIHI we do not subtract background 
but use 1.5% as the uncertainty from the background in 
the branching ratio measurement. Reweighting is used to 
correct the PID and momentum cut efficiencies. When 
the weight curve is changed by ler, the branching ratio 
changes by 0.3%. Different wire resolution models used 
in the MC produce shifts of 3.6% in the branching ratio. 
When the channel J/ip — > pp is simulated, the a value 
is an input parameter. If we change the value by ler, 
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FIG. 2: The (a) angular distribution for MC events; (b) total weight curve; (c) angular distribution of the background; and 
(d) angular distribution of the data. 



the efficiency changes by 0.6%. The total number of J/tjj 
events is determined to be (57.70 ± 2.72) x 10 6 [H]; we 
use 4.7% as the systematic error from this source. 

Combining the uncertainties from all sources in 
quadrature, the total effect on the branching ratio of 
J/tjj — > pp is about 6.2%, i.e. 

Br( Jj^ -> pp) = (2.26 ± 0.01 ± 0.14) x 10" 3 , 

where the first error is statistical and the second system- 
atic. 

Table Hill lists the systematic errors from all sources. 



VIII. SUMMARY 

The channel J/ip — > pp is studied with 57.7 x 10 6 
J/ip jl3j events. The branching ratio and angular dis- 
tribution are measured. The branching ratio is within 
one a of the PDG world average value and the angular 
distribution a is in good agreement with the theoreti- 
cal prediction that includes the effects of non-zero quark 
masses [8|. 



TABLE III: The systematic errors 



Sources 


Effect on a 


Effect on Br. 


MC statistics 




0.6% 


background 


2.2% 


1.5% 


wire resolution 


5.2% 


3.6% 


reweighting 


2.6% 


0.3% 


a value 




0.6% 


J/ip number 




4.7% 


total effect 


6.2% 


6.2% 
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